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Abstract: 



In order to elucidate the enhancement of pion production that may occur during dynam- 
ical scenarios of interest in connection with disoriented chiral condensates, we study the 
evolution of boson modes whose frequencies have a given arbitrary time dependence. The 
quantum-field treatment yields expressions for the time dependence of observables that de- 
pend only on the initial state and specific state-independent enhancement coefficients which 
can be obtained from the given evolution of the frequency. It is shown how these coefficients 
also can be obtained within an entirely classical framework by judicious usage of the depen- 
dence of the resulting enhancement on the phase of the initial state. Illustrative application 
is made for frequency evolutions obtained from idealized simulations of high-energy heavy- 
ion collisions and it is shown how the non-equilibrium evolution of the order parameter may 
lead to significant enhancements in the final pion spectra. 
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1 Introduction 



Several years ago, it was suggested that a high-energy colhsion may lead to a tran- 
sient approximate restoration of chiral symmetry within the very excited collision 
volume and, furthermore, the subsequent non-equilibrium relaxation of the chiral or- 
der parameter might lead to long-wavelength isospin-directed oscillations that were 
denoted as 'disoriented chiral condensates' (DCC) [|l], ^ |, ^, |^ . Effective field theories 
present an important theoretical tool for capturing the essential physics of this novel 
phenomenon, as the underlying quantum-field theory is too complicated to apply di- 
rectly. But even at the simpler level of an effective theory, the scenarios encountered 
in heavy-ion physics are very challenging because the dynamical evolution of the col- 
liding system requires that the field equations be solved in real time and, moreover, 
the system is often far away from equilibrium. 

The most commonly employed effective model is the linear a model . It can be 
treated in an approximate manner by means of the Hartree factorization technique, 
which linearizes the field equation and yields a description in terms of effectively 
independent quasiparticles subject to a common time-dependent effective field (which 
may be determined self-consistently from the local quasiparticle distribution). The 
evolution of the field, whether self-consistent or prescribed, then renders the effective 
quasiparticle mass time dependent. Thus, the physics of the quasiparticle modes 
corresponds to that of effectively free particles endowed with given time-dependent 
frequencies. Instructive developments of particular relevance to high-energy heavy-ion 
collisions have been made in recent years 0, ^, |^, |T0[ . 

The time dependence of the quasiparticle frequencies may generally enhance the 
pion yield. While the production of field quanta by a time-dependent field is well 
known in electromagnetism, the analogous phenomenon in harder to produce and 
explore in strong-interaction physics. An extreme example, the temporary creation 
of a supercritical field, may occur in a high-energy nucleus-nucleus collision if the 



hot interaction zone cools down sufficiently rapidly to generate a quench [11, 12 



If indeed such a quench occurs, then spontaneous pion production will be possible, 
in analogy with the electromagnetic scenario. But an accurate description requires 
quantum-field theory, as the vacuum fluctuations play an essential role as seeds of 
amplification. More generally, the production of particles by any time-dependent field 
cannot be accurately described within the realm of standard classical physics. 

This poses an important calculational obstacle, as the most tractable numerical 
simulation methods are basically classical in character and thus not adequate for 
an accurate calculation of the enhancement phenomenon. Thus, there is a need 
for further development of the dynamical treatment of the effective field theory, so 
that particle production by the time-dependent effective field can be described. The 
present study is intended to be a contribution towards this goal. As a step in this 
direction, we consider the simple situation when the spatial variation can be ignored, 
as may be justified in the bulk of the collision zone. With the restriction to idealized 
matter, the effective quasiparticle masses are independent of position and depend 
only on time and the discussion then simplifies significantly. 

The present paper has several distinct parts. First the needed quantum-field 
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formalism is introduced and the quasiparticle problem is brought onto diagonal form. 
Subsequently, a method originally conceived by Combescure [1^ is invoked to treat 



the evolution of each separate mode and it is shown how the development of various 
observables can be obtained simply in terms of specific state-independent coefficients. 
This result brings out the equal role played by vacuum and statistical fluctuations and 
elucidates how particle production can occur even when the initial mode is unagitated. 
Then we show how the same result for these coefficients can be obtained entirely 
within a classical treatment by means of a suitable decomposition of the dependence 
of the observables on the phase of the initial state. This result may prove useful 
for attempts to address the problem within classical field theory which is practically 
more tractable. Finally, the effect of the enhancement mechanism is illustrated for 
idealized dynamical scenarios of interest in connection with the DCC discussion. It is 
shown that typical time dependencies of the effective pion mass may lead to significant 
spectral enhancements. 

2 Framework 

In the Hartree approximation, the Hamiltonian density operator is given by 

7i:(x,t) = i[v^(x)t^(x) + V0(x)t.V0(x)+0(x)V(t)0(x)] . (1) 

The field operators 0(x) and '?/'(x) in (|I|) represent the local field strength and its 
time derivative, respectively. Since those are real, the associated field operators are 
Hermitean. Therefore, in three spatial dimensions, they can be Fourier expanded as 
follows, 

[flke"^-" + 4e-*'^-"] , (2) 

V^W = ^Ey^ 'a^e^^--aie~^^-] , (3) 

where fl is the volume. The operator aj^ creates a free particle having the wave vector 
k. Its frequency c^k is given by ujI = k"^ + rn?, where m is its free mass, and thus 
depends only on the magnitude of the wave vector, k = \k\. The free annihilation and 
creation operators satisfy the usual commutation relations for bosons, [a^^, aj^,] = 5kk', 
and it is elementary to verify the standard commutation relation for the field opera- 
tors, [0(x), '?/'(x')] = z5(x — x'). Finally, the total energy of the system is represented 
by the Hamiltonian operator if (t) = / d^x 7i(x, t). 

We may introduce the following "pseudo-canonical" variables, 

(4) 
(5) 
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They satisfy the usual canonical commutation relations, [Qic-Pk'] = ^^kk' and the 
total Hamiltonian can be rewritten as 



m = T.\[plK+^imiQv] , (6) 



where uj^ = k'^ + fi'^{t). However, the variables Qk and Pk are not hermitean, as is 
evident from eqs. (|^-|^), and they are therefore not fully suitable for our treatment. 

In order to bring the Hamiltonian on a truly diagonal form in terms of Hermitean 
canonical variables, it is necessary to perform a unitary transformation among any 
pair of modes having opposite momenta. In fact, there is an entire family of such 
transformations and we find it notationally convenient to adopt the following partic- 
ular one, 

f a^=k = TjK + a^k] 
\ aK=-k = -^[a-k-a^-k] ' 

where we have arbitrarily assigned specific opposite signs for the two members of each 
(degenerate) pair of opposite momentum labels k and — k. While the results will not 
depend on the particular convention, the explicit expressions for the transformation 
do depend on the assigned sign of the particular wave vector. In order to emphasize 
these new diagonal modes do not carry a directed momentum, they are labeled by 
the greek index k rather than k. Moreover, we will use = Uk to denote the free 
frequency of the two degenerate diagonal modes k = ±k. 

These new operators satisfy the usual commutation relations, = 5^^,', 

and the corresponding canonical operators are now Hermitean, 

= t , (8) 




= Pi , (9) 

with [qK^Vn'] = i^KK'- We also note the form of the operator representing the number 
of free quanta in the diagonal mode k, 



The effective Hamiltonian is now fully separated into independent harmonic modes. 
Hit) = EnH^it), where 

Hn{t) = \Pl + \u:l{t) ql , (11) 

so each mode k, can be treated separately. 

It is instructive to note that the field operators (0-0) can be rewritten as follows, 

^ [gKCOs(K ■ x) + g_^sin(K • x)] , (12) 

K,>0 



2\^ 

= ( — ) ^ [pkCos(k ■ x) +p_Ksin(K ■ x)] . (13) 

K>0 
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Consequently, the coordinate operators q^^ represent the coefficients in the trigono- 
metric expansion of the field strength 0(x), while the momentum operators repre- 
sent the coefficients in the corresponding expansion of its time derivative '?/'(x); those 
with positive momentum index go with the cosinus terms and the others do with 
the sinus terms (this particular feature is a consequence of the arbitrary enumera- 
tion convention introduced with the definition in (|^)). Of course one could make the 
trigonometric expansion in many other ways. This arbitrariness refiects the degener- 
acy of the modes with opposite values of k and corresponds to the freedom associated 
with a rotation in the space spanned by those two modes. 



3 Treatment 

After the diagonalization of H{t) has been acheived as described above, the evolution 
operator factorizes correspondingly, 

Uit.to) = \{U^{t,to) . (14) 

K 

It is then possible to treat each mode k, separately, using the method employed in ref. 



T3[| for a single mode. For each mode k, it is convenient to introduce the following 



two-particle operators, 

^0 = + > = , = . (15) 

They satisfy the following SU(1,1) algebra 



[Ko, K^] = ±K± , [K-,K+] = 2Ko , (16) 

and we note that Kq=Kq and K^_=K^. The Hamiltonian for the mode k, can then 
be reexpressed, 

H^{t) = 2uj+{t)k, + uj^{t)[k+ + k^] , (17) 

where uj±{t) = (ujlit) ±ujI)/2uJk for each k. 

It is now possible to define the following queeze operators that depend linearly on 
a time-dependent complex parameter T^it), 

A^it) ^ T^k+-Tlk^ , A^it) ^ T^k+ + Tlk_ . (18) 

We note that Al_=—A^ and A^^=Ak,- If we write F^; = Gf^exp^i^^^) we may instead 
use 7k = (y'K,exp(i^9^j;) as the parameter, where g^. = tanh(GK) (ensuring g^, < 1), and 
we have 

A = G[k+e'^ - k_e''^] , A = G[k+e'^ + i^.e"*'^] , (19) 

for each mode n. We may also introduce the squeezed annihilation operator (satisfy- 
ing = 

(3^ ^ e^-a^e-^- = [I - glV'^a^- ^^a^) . (20) 
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According to the method developed by Combescure [0, the dynamics of each 
of the decoupled modes k can be described in terms of a dimensionless complex 
parameter ^.nit) whose temporal evolution is governed by 

L + ujlit) = , (21) 

with suitable initial conditions (see later). Once ^K.{t) is known, the following auxiliary 
quantities follow, 

z^(t) = uj^^^ - , = aTg{z^) , (22) 
z^it) = u^^^ + i^^ , £^ = arg(z^) . (23) 

It should be noted that z^ is not equal to the complex conjugate z* . The parameter 
7k, is then given by 

7.(t) = ^^.e^'^'^ = ^ = ^ f , (24) 

SO that its phase is given by dn = ^k.~ ^k- The time evolution operator for the mode 
is then [|13| 

U^{t,to) = e^(^-W) e-2'("-(*)"^-(*"))^° e-^(^-(*o)) . (25) 

A number of simple cases having a constant frequency are especially instructive. 
With a view towards the later applications, we reserve the notation for the fre- 
quency of free motion, uj"^ = + m? , while tu^ generally denotes the frequency given 
by ujf^ = + yU^ where /i^(t) is the time-dependent effective mass squared. For nota- 
tional simplicity, we shall usually drop the mode index n and also often write ujq for 
the free frequency. 



3.1 Illustration: Free motion 

For free motion, /i^ = m? > 0, we have uj^^ = uf^ = + at all times. A simple 
solution to ( PI| ) is then ^nif) = exp{iujt), corresponding to 7K;(t) = and e^it) = u^t. 
In this case, A^, vanishes and the evolution operator becomes 17^, = exp{—2ieKo) = 
exp[—iH^t), which is recognized as the free evolution operator. Furthermore, the 
squeezed annihilation operator (pO| ) is identical to the free one, = a^. 

We also note that a coherent state of the form (p6|) , \x) = exp(xa"'' — x*a)|0), sat- 
isfies the free Schrodinger evolution equation, idt\t) = H\t) (with H = {a\a}u!o/2), 
iff the coefficient x(t) satisfies ix = ujqX (which implies x + ^oX = 0)- 



3.2 Illustration: Effective free motion 

At finite temperature, the application of the Hartree treatment to the linear a model 
leads to an approximate decoupling into modes that are described by an effective 
mass, /i^ > [0, H, H, [ly, |lT], |15|]. Then c<j^ = + /i^ is positive and time independent 
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and we may again seek a solution to (^) that leads to a constant 7^. Since we require 
g,. < 1 (recall = tanh G^) , we must choose 



7. = ^.e^''- = -"^^^ , (26) 

which is still real. But it is negative when Uk > So^, as is ordinarily the case in 
our applications (because the temperature increases the mass), and its phase is then 
'd^. = TT. For either sign, we have e^. = uj^t, corresponding to C,^ = exp^iu^t). The 
evolution operator then takes on the expected form, 

U^{t,to) = e-^'^^"-"-^* e-^ = e'^'^^'^-'^-^* , (27) 

analogous to that of a free particle with the mass fi, where the squeezed operators 
are those introduced in eq. (pO|). 

If the initial state is taken to be a coherent state of squeezed quanta (see Eq. (pGf ) 
in the Appendix), then the number of effective quanta is given by i> = {(5^(5) = Ixl"^, 
whereas the number of real quanta depends on the phase fj, u = (a^a) = |xP = 
(co'o/co')(Rex)^ + {uj/uJo)(lmx)'^. Moreover, the energy of the state is 

E = {H) = cu+(x*X + i) + + = {'y+'^)uJ+ + uu^cos2r] , (28) 



as follows directly from (p^Tf), and we see that E ranges between ulu (obtained if x is 
real) and uujq (obtained if x is imaginary). 



3.3 Illustration: Supercritical motion 

In this case we have fi^ < and we consider those (sufficiently soft) modes for which 
ft^ = —K^ — fi"^ is positive. We can achieve the simple initial value 70 = by adopting 
the solution = exp{Q^t) — 7+ exp(— fi^t), where 7+ = {luq + iVL,~)/{ijJo — iVL^). This 
leads to 

= ^+ ^2n^t_^2^ ' e^{t)-eo = aig[-^—^j . (29) 

Thus, for large times the modulus g^ approaches unity and, as we shall see in Sect. 
5.1, the particle number and related observables exhibit an exponential growth. 



4 Evolution of observables 

In the general case when fi^ evolves in time, we need to solve the equation (^) for 
each mode k, which is a straightforward task, once the initial conditions and /i^(t) are 
given. In scenarios encountered in high-energy nuclear collisions, the initial state is 
often described in terms of medium-modified quasi-particles with a suitable thermal 
effective mass. If the systems expands and cools, as is expected be the case, then /i^(t) 
decreases rapidly and usually approaches its free value in an oscillatory fashion. In the 
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most extreme cases, it may even attain negative values for a while. We wish to study 
such scenarios by employing forms of that have been obtained in dynamical 

simulations with the linear a model. 

Let us first consider the number of squeezed quanta in a given mode k, as counted 
by the effective number operator z/^ = [^]^13^, where we recall the definition of the 
squeezed annihilation operator (pO|), I3{t) = exp{A)aexp{—A). The expected number 
of squeezed quanta is then easy to calculate, 



Aog2jei^oQ,tQ, g-2jei^Og-Ao 



{to\PlPo\to) 



e e "e 



to) 



to) = (tole^oa^a e-^«|to) 



(30) 



where we have used the fact that the number operator a'^a commutes with Kq. 
The same calculation holds for any power (/?"''/?)". Thus the number distribution 
of squeezed quanta remains constant in time and therefore it provides no information 
on the dynamics. 

Using the relationship (|20| ) between a and P, we find 

1 



Therefore, if we start from a coherent state of the form described in the Appendix 



(31) 



where x^,= \Xk\ exp(-2r7^) 



we have 



2,1 
^ 2 ' 



2(to|i^±|to) 



Consequently, the initial number of physical quanta is 

fQ = (tol^'lto) = (i^o|2-ft'o|^o) = 



while, with 70 = (70 exp(i'i9o)) the initial number of squeezed quanta is 

1 



z/(to 



1 



9l 



1 + ^0 - 2 cos(^9o + 2r7) z/o + ^0 



(32) 
(33) 
(34) 

(35) 



To extract interesting dynamical information, it is useful to note that the subspace 
spanned by the three rotation operators Km is closed under the time evolution, 

kmit) = U^{t,to)KmU{t,to) = J2Zmrn'{t) Km' . (36) 



The complex time-dependent structure coefficients Z„ 



mentary means. Writing Z„ 



can be determined by ele- 
we find 



^00 



:i + (7o)(l + ^?') - 4(70(7 cos A , 
-(7o(l + 9^) + (1 + 9o)9 cos A + z(l - gl)g sin A 



-2(70(7 + (l + (7or) cos A + z(l 



^^/)sinA 



+ 9o)9 - 9oi^ + 9 )cosA-z5(o(l -g )sinA 
-2gog + {gl + g'^) cos A - i{gl - g^) sin A 



e^''« = 

j(f-i5o) 



0-) 



z 



0' 



^j(i9+i?o) 



z 



+-' 



(37) 
(38) 

(39) 

(40) 

(41) 
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where A = e + e — Eq — Eq and we recall 7 = gexp{{}). We note that the diagonal 
coefficients start out from unity, Zmmito) = 1, while the off-diagonal coefficients 
initially vanish. Moreover, Zqq can never become smaller than one, 

{l + gl)il + g^)-4gog _ {1 - g,g)^ + {g - g^f 

The evolution of any observable quantities should not depend on the choice of initial 
values ^0 and ,^0) because any choice leads to the same evolution operator U. This 
physical feature is reflected in the fact that the above Z coefficients are in fact all 
independent of the initial conditions ^0 and ^o- Furthermore, it is elementary to show 
that Zoo(O) = uj^, which shows explicitly that the initial growth of Zqq is independent 
of ^0 and ^0 and depends only on the physical quantity uj_. We flnally note that the 
magnitude of each Z coefficient remains constant during time intervals where u'^ is 
unchanged. 

That the fact that the functions Zmm'{t) are independent of the specific initial 
values of ^0 and ^0 implies that they depend only on the magnitude k = |k;|, since 
the equation (|2T| ) for contains only cu^. Thus, in effect, Zmm'if) depends only on 



the initial energy ujq. 



4.1 Expected particle number 



We now focus on the development of the expected number of free quanta in each of 
the independent modes, i^K,(t) = (t|a^a^|i(:) (keeping in mind that this quantity is 
distinct from the number of quanta having a given momentum, ^^(t) = {t\a\xi-^^\t)) . 
Since the number operator can be expressed as z> = = 2Kq — |, we find 



u{t) = U^uU = 2Zoo{t)ko + 2Zo+{t)k+ + 2Zo_(t)i^_ 



(43) 



If we start from a coherent state of the standard form (|86D , with x=yi^ exp(— zr/), 
we easily find 



Zoomi^o + l) - I + [^o+(t)e^** + Zo-(t) e-— jz/o 



-2ir?ol 



9l) 



- g'r' {[(1 + gl)il + - Agog cos A]( 



^0 



(44) 



- 2 



go{l + g"^) cos(2r7o+?9o) - 9 cos(2?7o+?9o+A) - g^g cos(2?7o+?9o-A) z^o} 



We note that if we start from vacuum then z/q = and the resulting number of quanta, 
z/(t), is given solely in terms of the coefficient Zqq, v = \{Zqq — 1), which does not 
contain the initial phase 779. This is in accordance with our expectations, since the 
vacuum has no definite phase. In the present description this feature is automatically 
guaranteed by the fact that the coherent-state parameter for the vacuum vanishes, 
Xvac = 0, and so the phase has no significance. 

If we consider an ensemble of coherent states for which the initial phase rjQ is ran- 
dom, -< e*''" )^=0, then the ensemble average of ( PD eliminates the terms containing 
r]Q. [We note that this phase average can be obtained exactly by averaging over N 
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equally spaced angles rjn = riQ + 27m/N, n = 1, . . . ,N, where N>2 and the value of tiq 
is immaterial.] The evolution of the particle number is then again given by the first 
term in expression (0), which may be recast to exhibit the increase in the occupancy, 



-< Au{t) y = -i v{t) y -uo = (Zoo(t) - 1) {uo + |) • (45) 

That expression has a simple interpretation: the field fluctuations, which are initially 
proportional to z/q + |, are magnified by the factor Zoo(t). The total field fluctua- 
tions are represented by the operator 2Kq and contain both the vacuum fluctuations 
(amounting to half a quantum for each mode k) and the fluctuations associated with 
the real excitations of the mode (represented by the number operator 0). The above 
result may also be written on the form 

-< u{t) + 1 ^ = XN{t) ^ u{to) + I ^ , (46) 

where we have introduced the number amplification coefficient Xn = Zqq. It follows 
from our earlier discussion that this amplifcation coefficient depends only on the mag- 
nitude of the wave number, which in turn implies that the original oppositely moving 
states ±k are on the average affected in the same manner by the time dependence 
of w^. However, this is not true for the individual states due to the phase-dependent 
terms in (|4^). We also note that the above expression corresponds to the result 
derived by Boyanovsky et al. in the Hartree approximation [|, |T0[] . 

4.2 Variance in particle number 

We turn now to the variance in the number of quanta in a given mode, 

^l{t) = {t\ul\t) - . (47) 

Because of the subtraction, we need only consider the commutator terms resulting 
from bringing the square of the number operator on normal form. We therefore find 

a\t) = (to|^oV + 4^o(^o+i^+ + ^o-i^-)+4Zo+Zo_(z>+i)|to) (48) 
= + 2ZooiZo+e^'^ + Zo-e-^'") + 4Zo+Zo_] uo + 2Zo+Z^^ . (49) 

If again we perform an ensemble average, the terms with the phases t] disappear, 

^a\t)y= (Z2o + 4Zo+Zo_)z/o + 2Zo+Zo- = X^z/q + 4|Zo+|2(z/o + |) • (50) 

Since .^oo(^o) = 1 aiid 2'o±(to) = 0, we initially have cr^(to) = ^(to), as is characteristic 
of the Poisson distribution implied by the coherent state. But, except when the 
motion is free, the value of o"^ will generally differ from u for t > to, and hence the 



time-dependent state does noi retain a standard coherent form (c/. Ref. [[14| ). 
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4.3 Expected energy 

Finally, we consider the evolution of the energy expectation value in a given mode, 

E{t) = {t\H\t) = {t\2uo+Ko + uo-{k+ + k^)\t) , (51) 

where we recall that uj± = (cu^ — ujI)/2ujq. Again, we may use the relation ( p6D to 
express the time- dependent expectation values in terms of expectation values in the 
initial state, 



E{t) 



[uj+Zqq + u;_Re(Z+o)] («^o 



(52) 



+ 



2c^+Re(Zo+ e^^'"') + cu.Re ((Z++ + Z_+) e 



Furthermore, for an ensemble in which the expectation values (to|-^±|^o) = t'o exp(±2z?7o) 
average to zero, only the first line of (p^) contributes and we obtain 

^ E{t) y = {u+{t)Zoo{t) + c^-(t)Re(Z+o(t))} {i^o + k) = Xsit) ^ E{to) y 



/)-! {uj+ [(1 + gl){l + g^) - 4(70(7 cos A 



2cu- 



(53) 

;i + gl)g cos d-go co^{d - A) - g^g"" cos{d + A)] } (z/q + \) . 



where Xe = [uj+Zoo+Lj^Re{Z^o)]/uj+{to) is the energy amplification coefficient. With 
e = E/u^ito) — |, the above result may be recast in a form analogous to (^6|), 



-< e{t) + i ^ = XEit) -< e{to) + i ^ 



(54) 



We note that when the frequency has returned to its free value, u"^ = Uq (as usually 
happens sooner or later), then U- = and uj+ = uo- Thus the expression (|5^) for 



the energy enhancement becomes identical to expression (|4^ for the enhancement of 
z/ + |, as one would expect since then E = + ^)uo. 

For motion with a constant positive value of u, the quantity 7 remains constant 
in time, 7 = 7_ = — (tu — u!o)/{u! + uq), and so we have 



XN{t) 



Zm{t) 
(t) + Z_+(t) 



^--^cos2ujt > 

.UJ- 



1 , 



1 — COS 2ujt] + i — sin 2ujt 

,2 



-2Z 



0± 



COS 2ujt 



I — sm2ci;t . 



(55) 
(56) 

(57) 



Then the time-dependent terms in (|5^) cancel out, as one would expect, but the 
resulting conserved energy depends on the initial phase tiq ii ui ujq, 

E{t) = E{to) = {uo + \)lo+ + z/o o;_cos2r/o , (58) 

and obviously the energy amplification coefficient remains equal to unity, Xs(t) = 1. 



10 



5 Classical treatment 



We have described above how the cranked harmonic oscillator can be treated exactly 
with the tools of quantum field theory and we have illustrated how the development of 
observables can be obtained by means of appropriate state-independent amplification 
coefficients whose evolution depends only on the prescribed time dependence of the 
frequency of the mode. We now describe how the key results can be obtained alterna- 
tively within the framework of classical dynamics, since this may be practically useful 
in the context of dynamical simulations. We shall restrict the present discussion to 
scenarios where the initial frequency is the free one, u!^{t < 0) = uJq. 

The starting point is the expression ( p^T]) for the Hamilton operator H in terms of 
the canonical operators q and p which suggests considering the classical Hamiltonian 
function 

n{q,p) = + \uj\t) . (59) 

The corresponding equations of motion are the familiar Hamilton equations, 

q = p , p = —uj'^{t)q . (60) 

For a given initial quantum many-body state of standard coherent form, |xo); we 
associate a specific classical phase point {qo,po), where 



\J (2z/o + l)/ujQ cos?7o , = -\J (2/^0 + l)t^o sin?7o . (61) 



We recall that the complex coherent-state parameter is given on the form xo = 
y^exp(— ir^o), with = |xoP being the expected number of quanta in the coherent 
state, z/q = (xo|ct^Q^|Xo)) and ujq being the initial (free) frequency which is positive. We 
then solve the classical equations of motion (|60D, using {qo,Po) as the initial condition 
and with the same specified uj'^{t) as employed in the quantal calculation. 

The time-dependent energy of the classical system, S, is then given by the value 
of the Hamiltonian function, 

S{t) = n{q{t),p{t)) = ip(t)2 + y{t)q{tf . (62) 

Furthermore, in order to extract number of quanta, we note that 



and we therefore define the corresponding classical observable, 

m ^ ^ pitf + ^ qity . (64) 

It is seen from eqs. ( ^ ) and (|5^ ) that the quantal result for the time-dependent 
expectation values of particle number and energy, u and E, depend in a very sim- 
ple manner on the initial phase tjq, each being a constant plus a second harmonic. 
Therefore we may write 

u{t) = -< u{t) >- + 6u{t) cos(277o - Sn) , (65) 
E{t) = -< E{t) >- + 6E{t) cos(2r/o - 5e) ■ (66) 
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Since the basic equation of motion is the same, when expressed in the canonical 
variables, we expect a similar behavior of the classical result. This is indeed borne 
out by our calculations and we may therefore write the classical results on a similar 
form, 

X{t) = -< Ar{t) >- + SAfit) cos(2r/o - S'j^) , (67) 
n{t) = -< H{t) >- + 6n{t) cos(2r7o - S'^) , (68) 

where the phases 5^ and 6'^ governing the classical results are found to be the same 
as those in the expressions (|65|-p6|) for the quantal case. 



In the quantal case, the coefficients in (|65| - pq ) are given in terms of the calculated 
Z coefficients. In the classical case, the coefficients in ( ^7| - p8| ) can be extracted by 
Fourier resolving an ensemble of dynamical results with respect to the initial phase 
rjo. Since only the zeroth and second harmonics are present, it suffices to consider 
three separate dynamical histories having their initial phases separated by 27r/3. This 
elementary analysis then yields the evolution of the three parameters in each of the 
two expressions ( |67D and (|68|). 



The quantal result (p5|-|5BD may also be written as 



u{t) + 1 = (z/o + l)Xjv(t) + jyoSX^it) cos(2r7o - 6^) , (69) 
E{t)/uo = iuo + ^)XEit) + iyo5XEit)cosi2r]o-6E) . (70) 

It is then seen that -< u >~= (z/q + j)Xn — ^ and -< E >~= (z/q + ^)u!oXe, as we 
have found already. Furthermore, 6^ = u^dX^ and 5E = uqUoSXe- Thus the phase- 
independent terms are proportional to i^o + |? whereas the phase-dependent terms are 
proportional to z/q, which is a reflection of the presence of vacuum fluctuations in the 
quantal ground state. 

By contrast, in the classical treatment all terms scale in the same manner, being 
proportional to z/q + A particular classical calculation is therefore only a quanti- 
tatively good approximation to the quantal result when the difference between z^o + | 
and z/q is immaterial, i.e. for large occupation numbers, z/q ^ 1. However, as the 
above analysis shows, the described resolution of the phase dependence has enabled 
us to separate the two terms in the classical treatment and we can therefore readily 
define the corresponding amplification coefficients, 

_ ^Afy , _ 5M 

= — ) = — -T ' y''-) 

= 7 — ' ^-^^v = 7 -T\ — • C^^) 

(z/q + ^)UJo (z/q + 

Once these coefficients have been determined on the basis of the classical calculation, 
by analysis of three separate trajectories, they can then be used in the expressions 



^The specific choice of i^o + 5 for the action in the classical calculation is not essential due to the 
universal scaling; the only important thing is to always employ a finite value for the action, so that 
even the evolution of the quantal vacuum can be addressed. 
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( |69D and ([7^) to provide the evolution of the quantum expectation values z/(t) and 
E{t) for any initial phase tiq and any value of the initial number of quanta uq. We 
emphasize that this method provides the exact quantum result even though it has 
been obtained entirely by classical means. 

The above method provides an alternative means for obtaining the values of the 
amplification coefficients X^v and Xe solely on the basis of classical mechanics. In 
order to demonstrate the utility of the method, we now consider a number of especially 
simple cases. 



5.1 Test case: Constant negative cj^ 

We first consider the case when a free mode is suddenly exposed to a supercritical 
field. So V, 
easily find 



field. So we assume u;^ = cUq > for t < and cu^ = —Q"^ < for t > and then 



Xiv(t>0) = 1 + 2 f "^"^!^' ^ sinh^fit. (73) 

Thus if we start from vacuum, z/q = 0, the number of spontaneously created quanta 
grows approximately exponentially while the supercritical value u"^ = —Q"^ is main- 
tained, 

u, = 0: u{t>0) = |(X^ - 1) = f si^h^ ~ e'""* • (74) 

To obtain the number amplification coefficient with the corresponding classical 
treatment, we use the form q = aexp(f2t) + bexp{—Qt) and readily find the initial 
phase point. 



qo = a + b= J— cosr/o , Po = ^{a - b) = - J (2i^o + l)t^o sinr^o , (75) 

V ^0 

from which the expansion coefficients are readily determined. The time evolution is 
then 



q(t)= M-[e^" + e-^"]cosr7o + ^-^[e^« - e"^"] sin r/o| , (76) 



(z/Q + |)fi f r^rm f^i^nt , -nti 



p{t) = f 2 \V^[' ' ~ \I^K' + ^-n^^^Voj . (77) 

The time evolution of A/" = p"^ /2uo + q'^Uo/2 then follows and we note that A/q = ^o + ^- 
If we average over the initial phase t/q the mixed terms vanish, leaving 



XL = r- = r -< P + —Q >~ = 1 + 2 — — — 



2 I n2\2 

sinh^ Qt , (78) 



which is the exact same result ( |73D as obtained above with the quantal treatment. 
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Figure 1: Time dependence in supercritical scenario. 



The time dependence of the number of quanta in the system, i^{t), for a simple supercritical 
scenario in which uj{t < 0) = loq and uj'^^t > 0) = —uJq. The number of free quanta as a 
function of the elapsed time t since the sudden onset of supercriticality, averaged over the 
initial phase ijq. The values of the initial number of quanta vq are indicated. Solid curves: 
quantum calculation; solid dots: classical calculation. Corresponding quantal and classical 
results are identical to within the numerical tolerance. 

The resulting growth of the number of quanta is illustrated in fig. || for the par- 
ticular case where Q = ujq. The results have been obtained by numerical solution of 
the respective quantal or classical equations of motion, as described earlier, and the 
two treatments give the same result, = X^r, as the above derivation demands. 

For the evolution of individual states, the agreement between the quantal and 
classical calculations is illustrated in fig. | which displays the number of quanta 
present after tu^ has maintained its constant negative value — cJq for a certain length 
of time. At = 1 fm/c. Since the phase rjo is a cyclic variable, the results are periodic 
in this parameter, but the period is only vr since the phase dependence enters via the 
combination 2?7o. It should be noted that the two treatments are in full agreement 
as well with regard to the phase- dependent results, thus demonstrating that also 
the quantities governing the phase dependence are identical for both amplitude and 
phase, SX'j^ = 5X^ and 5'^ = 6]y. 

5.2 Test case: Constant positive cj^ 

When u = uo for t < and u = const > for t > 0, the quantal result is 



XN{t > 0) 



^ + — cos2cut = 1 + 




(79) 
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Figure 2: Phase dependence in supercritical scenario. 
The number of free quanta v{tj) at the time = 1, as a function of the initial phase 
r/o, for the same supercritical scenario considered in fig. ||. The initial number of quanta 
vq is indicated. Solid curves: quantum calculation; solid dots: classical calculation. The 
corresponding phase-averaged result is shown by the horizontal dashed line for each non- 
zero value of vq. Corresponding quantal and classical results are identical to within the 
numerical tolerance. 



To derive the classical result, we use the form q = a cos tut + b sin cut and readily find 
the initial phase point, 



a = \I2vq + 1/ujq cosrjQ , po = ub = —\/1vq -|- Xloq sinr/o • (80) 



It is then elementary to obtain the amplification coefficient, 

^^ ,ft_if ^= cos2^t+l(^4 + 4)si^'^^ = (81) 

Thus, also in this case, the quantal and classical expressions are identical. 

We note that the number amplification coefficient oscillates between its initial 
value Xfqit = 0) = 1 and its maximum value X^^'^ = {u"^ + uJq)/2lo'^Uq, and the 
temporal average is given by X^"^ = [(tu^ + ci;Q)/2ti;ti;o]^- The amplitude increases with 
the difference between u and Uq and the frequency of the oscillation is equal to 2lj. 
Therefore, the resulting number of quanta in the system depends sensitively on the 
duration of the perturbation. 

5.3 Test case: Periodic positive cj^ 

Finally, we consider the special case of a harmonic perturbation of free motion, 

uj'^{t) = uj^ [l + AsinQt] . (82) 
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The equation of motion then has the Mathieu form and so we expect a resonant 
behavior as a function of the perturbation frequency Q. This is illustrated in fig. |^ 
for the number enhancement coefficient Xj\f resulting from subjecting the system to 
two periods of the harmonic perturbation. For small amplitudes and many periods 
these are given hj Q = 2ujQ/n, where = 1, 2, 3, ... is a natural number. The reso- 
nance frequencies shift downwards when the perturbation amplitude A is increased. 
Moreover, it is of course not possible to accumulate much amplification during only 
two periods unless the amplitude is large. While the low frequency region, fl < uq, is 
characterized by rapid oscillations, with an average X^v value of around eight when 
A is unity, the higher region exhibits one broad bump centered somewhat below 2ujq. 
This latter feature implies that a given perturbation frequency Q is especially effective 
in enhancing those modes that have about half that frequency, cj^ ~ It should 
be noted that the induced enhancement obtained with the classical treatment agrees 
with the quantal results. 
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Figure 3: Number enhancement in periodic scenario. 
The final number enhancement coefficient X^- after the free osciUator has been subjected 
to two periods of a harmonic perturbation, = cjq [1 + ^ sin Qt] , as a function of the per- 
turbation frequency 0, in units of wq; three values of the amplitude A have been considered, 
A = 0.50,0.75, 1.00. The curves have been obtained by means of the quantum calculation, 
X^^, while the solid dots result from the corresponding classical calculation, X^^. 

The dependence on the initial phase tjq is shown in fig. |^ for the large amplitude, 
A = 1, in the case when one initial quantum is present, z/q = 1. The behavior 
is qualitatively similar to what was shown in fig. ^ with the classical calculation 
tracking the quantal behavior perfectly. The average value of the final occupancy Uf 
is about 16, in accordance with fig. ^ where the amplification factor X{Q = 2ujq) is 
seen to be about 11 (indeed, Uf = (z/q -|- |)X ~|~|"11~| = 16). 
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Figure 4: Phase dependence in periodic scenario. 
The final number of free quanta Vf asa function of the initial phase r/o , after exposing the free 
oscillator to two periods of a harmonic perturbation, = Li;g[l + sin2u;o^] (corresponding 
to an amplitude A=l and a frequency 17 = 2ujo); the results are shown for unit initial 
occupancy, fo=l- Solid curves: quantum calculation; solid dots: classical calculation. The 
phase-averaged results are shown by the dashed horizontal line for each fQ (it is the same 
for Qu and CI to within the line width) . 

5.4 Semi-classical treatment 

In many numerical studies the field operator is replaced by a classical field, 0(x) —>■ 
0(x), and the equation of motion can then be solved by standard means (see refs. 

[1^ for examples involving the linear a model). In this type of approach, the 
field strength is taken to represent the degree of agitation, without allowing for any 
contribution from the vacuum. Thus the number of quanta present in a given mode 
K is taken to be proportional to the square of the corresponding field amplitude, 

~ This corresponds to using z/q rather than i/q + | as the action in the 

classical phase-space dynamics (see eq. (^1])). Thus, the quantum vacuum would be 
mapped into the ground state of the classical oscillator, |0) ^ {q = 0,p = 0), where 
it would then remain for any given form of c<j^(t), thus precluding the spontaneous 
creation of particles from the vacuum. 

Moreover, the final phase-averaged action would become = t'o^Af which is 
always smaller than the quantal result, z/^" = (z/q + ^)X]y — |, since Xj^ > 1- While 
the two results are practically similar for large actions, uq ^ 1, the typical situations 
encountered involve occupation numbers well below unity. In such cases the semi- 
classical treatment then leads to a significant underestimate of the degree of particle 
creation induced by the imposed external field. We shall return to this important 
point in the discussion of the quench scenario in Sect. 6.3 and a specific example is 
shown in fig. ^ 

The above described treatment makes it possible to circumvent these inherent 
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shortcomings of the classical treatment and extract results from the classical calcula- 
tion that are quantitatively accurate for any degree of occupation. 



6 Application to DCC dynamics 

We here consider forms of uj'^it) are especially relevant to scenarios of interest in 
the context of disoriented chiral condensates. They correspond to the conditions 
prevailing in the interior of initially hot but rapidly expanding matter. 

6.1 Initialization 

In the discussion of disoriented chiral condensates, it is usually assumed that the early 
part of the nuclear collision generates a very excited volume within which chiral sym- 
metry is at least partially restored. The isospin degrees of freedom are then described 
in terms of effective quasipions characterized by a medium-modified dispersion rela- 
tion, a;^ = k'^ + n'^{t). The effective mass depends both on the magnitude of the chiral 
order parameter and the degree of agitation which both evolve dynamically. The 
initial degree of agitation can conveniently be characterized by an equivalent temper- 
ature To which may significantly exceed the free pion mass rriT^. The corresponding 
effective mass /iq is then also significantly larger than m^. In our applications, we 
shall employ estimates obtained with the linear a model which has been the most 
popular tool for DCC studies. 

After its preparation in a state of high excitation, the system is expected to 
expand rapidly and the corresponding sudden drop in the effective mass may lead 
to quench in which the effective field has become supercritical, /i^ < 0. In order 
to generate such a quench, it is necessary to first achieve approximate restoration of 
chiral symmetry and this in turn requires a fairly high initial agitation, corresponding 



to an equivalent temperature of several hundred MeV [|T2|. To illustrate this we 
shall employ Tq = 400 MeV which ensures that the initial order parameter is small, 
(To < 10 MeV; the corresponding effective pion mass is taken as fio = 560 MeV/c^. 
(We recall that in vacuum the order parameter is equal to the pion decay constant, 
(To = /tt ~ 92 MeV.) Alternatively, in order to explore the possibility that such high 
initial excitations are in fact not attained, we shall also consider a more moderate 
value. To = 200 MeV, for which the order parameter is still fairly large, a ~ 60 MeV; 
the corresponding effective pion mass is taken as no = 240 MeV/c^. The specific 
parameter values are based on the semi-classical studies in ref. [|15| which lead to 
results that do not differ essentially to those obtained by other methods (such as the 
recently developed optimized perturbation theory |T8|). 



In fig. ^ we illustrate the relationship between the various occupation numbers 



of interest. We assume that the initial state is of the coherent form ( pUj) defined in 
terms of the effective (or squeezed) quasipion states. In the present discussion, we are 
interested only in quantities averaged over the initial thermal ensemble and then the 
results depend only on the magnitude of the state index, k = \k,\ = |k| = k, which we 
shall therefore replace by the magnitude of the pion momentum, p. (The momentum 
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Figure 5: Initial number of quanta. 
The solid curve shows the initial number of effective quanta vq in modes for which the 
momentum has the magnitude p, under the assumption of thermal equilibrium at the tem- 
perature To and using a corresponding effective pion mass [1q as given by the linear a 
model 1 15]; the dispersion relation for initial effective pion mode is then = + /ig- The 
dashed curve shows the corresponding number of real quanta, = VQiOA^juJi. The thermal 
occupation number for a free pion gas (//q = "i.^-) is also shown (dotted curve). 



is a constant of motion in the present scenario where all the dynamics arises from 
the time dependence of /i^.) The initial energy uji of a given mode is then given by 
the dispersion relation uJl = + l^'i- When expressed in terms of the quasipions, the 
initial Hamiltonian operator is given by Hito) = ^uJi{P\ P}, as we have seen, and 
the thermal average of the occupation number is therefore of the usual Bose- Einstein 
form, 

z>o = ^ iXolP^PlXo) y= [exp(^) - 1]-^ . (83) 

This quantity is always well below unity (in fact it never exceeds about one half at 
any temperature |]15[) and it decreases towards a constant value at high temperatures 
because the effective pion mass then grows in proportion to Tq). 

The number of real pion quanta is represented by the number operator a^a and 
it follows from our earlier analysis that we have 

i^o = -< (Xo|""^"IXo) ^ = — J^o >vo, (84) 

since 0;+ = \{ujf + uJq)/ujq. This quantity is always larger than z/q, especially for 
the softest modes where the modification of the mass has the greatest influence. It 
follows that the form of z/o(p) cannot be approximated by a free pion gas at a suitably 
adjusted effective temperature temperature. Finally, it is important to realize that the 
corresponding number for a non-interacting gas of free pions at the given temperature 
To would be significantly higher, due to the (unjustified) use of the free mass m-,,. 
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6.2 Evolution of the mass 



In order to emulate the conditions that may be experienced inside the expanding 
coUision zone, we employ two different forms for /^^(t) corresponding to two opposite 
idealized expansion scenarios. 

For the quench scenario we employ the time-dependent effective pion mass ob- 
tained in ref. by subjecting matter in initial equilibrium at Tq = 400 MeV to a 
Rayleigh cooling with D=3. (which emulates the effect of a three-dimensional scal- 
ing expansion 0). The resulting /i^ is shown in fig. ^. Starting with a high 
value (reflecting the inital thermal mass), it rapidly drops into the negative regime, 
followed by large but decreasing oscillations (the amplitudes drop off as ~ For 
convenience, the oscillations have been artificially suppressed above t ^ 6 fm/c. 




012345678 4 8 12 16 20 



Elapsed time f(fm/c) 
Figure 6: Initial number of quanta. 
The time dependence of the square of the effective mass, /^^(i), in the two scenarios 
considered. Left panel: The quench scenario obtained by first preparing the system at 
To = 400 MeV and then cooling with a rate emulating a three-dimensional scaHng ex- 
pansion; taken from ref. with the oscillations suppressed beyond t = 6 fm/c. Right 
panel: The moderate scenario in which five initial conditions are sampled from the thermal 
ensemble prepared at Tq = 200 MeV and then evolved up to t = 20 fm/c. 

In the more moderate scenario, where the systems starts from Tq = 200 MeV, we 
have assumed that the system experiences a standard longitudinal scaling expansion 
and solved the corresponding field equation by switching from (t, z) to (r, 77) (see, for 
example, refs. p, |I4[ |1^). Since the order parameter is initially quite large (about 
60 MeV), the system never enters the unstable regime and the effective mass never 
drops much below the free value. But the oscillations nevertheless persist for a long 
time (their amplitude drops off as ~ t~^/^) and we have truncated the evolution 
at t = 20 fm/c. In order to illustrate to effect of event-by-event fluctuations, we 
show evolutions that start from five different initial conditions, each one sampled 
from the thermal ensemble.^ Since the number of quanta in each individual mode 

^The interested reader is referred to ref. [|l5j for a discussion of how the initial state can be 
sampled from a thermal ensemble of semi-classical fields. 
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is stochastic (governed by a Poisson distribution), the resulting effective mass also 
exhibits ffuctuations. Nevertheless, the different evolutions lead to slowly subsiding 
and fairly regular oscillations which have approximately the same frequency. 

With the system initialized as described above, and with the time dependence of 
given, it is now straightforward to solve the equations of motion as described 
earlier and obtain the number enhancement coefficient X^v for each particular pion 
mode. We will discuss the two scenarios in turn. 



6.3 Quench 

The quenched evolution can be regarded as a combination of the supercritical and 
oscillating scenarios considered in sects. 5.1 and 5.3, respectively. The excursions 
into the supercritical regime lead to a significant degree of pion production. This 
effect is largest for the softest modes, since those have the smallest uP' . As a result, 
the amplification coefficient X^r grows steadily as p is decreased towards zero. The 
oscillations in [j? lead to the prominent peak near p = 200 MeV/c^, reminiscent of 
the broad peak in fig. ^. If the oscillations had not been truncated, this component 
would have become even more dominant. The resulting amplification coefficient is 
shown in Fig. |^ (left panel). It exceeds unity by about one order of magnitude up to 
p ?a 260 MeV, at which point it quickly drops to unity. 

The final yield of pions can now readily be obtained. For each value of p we know 
the initial number of real quanta vq (given by the dashed curve in fig. §), and the final 
number of free pions is then given by = (z/q + ^)Xiy — |, as discussed earlier. The 
resulting spectral profile is shown in the right panel of fig. ^. The initial values are 
about one half throughout the region affected (see fig. i.e. the thermal fluctuations 
are comparable to the vacuum fluctuations. The resulting enhancement of the yield 
is therefore roughly equal to twice the value of Xjv, which amounts to a factor of 
20-30 up to 260 MeV/c, in the particular case presented. 

In a semi-classical treatment of the field theory, the initial field amplitude does 
not contain contributions from vacuum fluctuations. Consequently, the final yield 
would be given by 



- uoX^ < (z/o + i)Xjv -I = uf\ (85) 



as already discussed in Sect. 5.4. This result is included in fig. |^. If the initial occupa- 
tion is large, uq ^ 1, then the semi-classical result is a good approximation. However, 
in the cases of present interest, the occupancies are usually at most about one half, so 
the semi-classical calculation should be expected to underestimate the enhancement 



by at least a factor of two. Thus, for example, the enhancements obtained in ref. |T2 
with the semi-classical treatment of the linear a model should be regarded as lower 
bounds on the expected effect. 

Finally, we wish to express the opinion that the present quench scenario is probably 
rather optimistic. In order for the quench to occur, approximate chiral symmetry must 
first be generated throughout the volume considered; this is usually thought to occur 
as a result of the expected high degree of agitation produced in the early part of 
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Figure 7: Result of the quench scenario. 
Left panel: The final value of the number amplification coefficient as a function of 
the magnitude of the pion momentum, calculated with the time-dependent mass associated 
with the quench scenario (left part of fig. Right panel: The corresponding final invariant 
differential pion yield. Solid curve: the final yield of free pions; dashed curve: the number 
of real quanta in the initial state; dotted curve: the final yield obtained with a semi-classical 
treatment. 



the high-energy collision. Subsequently, the system must cool sufficiently rapidly so 
as to recover approximately the low-temperature effective potential before the order 
parameter has relaxed back to its vacuum value. This delicate race requires a very 
rapid expansion. A variety of calculations suggest that the usual longitudinal scaling 
expansion is insufficient |jl2|, |1^ and it may be questionable whether the onset of the 
ultimate transverse expansion could occur before the inherent longitudinal expansion 
had already caused the order parameter to grow too large for a supercritical condition 
to arise. It therefore appears reasonable to assume that the quench scenario represents 
an extreme limit. 



6.4 Fizzle 

We now turn to the (probably more likely) scenario, in which the early dynamics 
does not achieve approximate restoration of chiral symmetry. As a specific example, 
we assume that the initial condition corresponds to a temperature of "only" Tq = 
200 MeV. Since this is still a bit below the range at which the order parameter 
decreases most rapidly (see for example refs. [|r^, [T^ for specific recent calculations 
within the linear a model), the equilibrium value of the order parameter is still fairly 
high, (To ~ 60 MeV, and even a three-dimensional scaling expansion would not be 



able to bring the system into the supercritical region Therefore, the primary 

enhancement mechanism is the relatively moderate oscillations in the neighborhood 
of the ground state. 

In the left panel of fig. H is shown the resulting values of the amplification coef- 
ficient Xn as a function of the momentum p, for each of the five histories based on 
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the initial conditions that were sampled from the thermal ensemble of quasiparticles 
at the given Tq. Because of the differences in the microscopic composition of the five 
initial states, the detailed evolutions will differ somewhat and, as a consequence, we 
will obtain an entire ensemble of resulting Xisi{p). However, although they differ in 
detail, the various histories tend to yield qualitatively similar results, as is evident 
from the display. 

The result can be characterized as a broad and modest increase of X^- by about 
20% up to about p ^ 240 MeV/c plus a relatively narrow peak centered at about 
that momentum, followed by an abrupt drop-off at higher momenta. The peak is 
positioned where we would expect: It follows from the discussion in Sect. 5.3 that 
the modes most affected by the quasi-regular oscillations of the mass are those with a 
frequency just below half the a mass (the calculations use m„ = 600 MeV/c^). Since 
the effective pion mass is about /iave ~ 150 MeV through the later part of the evo- 
lution, the frequency associated with the most affected mode is about u ~ 285 MeV 
which is indeed slightly less than ^rria- We note that in the present calculation the 
oscillations in the effective mass were artificially eliminated after tmax = 20 fm/c, 
even though the system would keep ringing for a longer time. Naturally, if tmax were 
increased then the resulting amplification coefficients would be larger. 
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Figure 8: Result of the fizzle scenario. 
Left panel: The final vaiue of the number amplification coefficient X]\f as a function of 
the magnitude of the pion momentum, calculated with the time-dependent mass associated 
with the fizzle scenario (right part of fig. |6|). Right panel: The corresponding final invariant 
differential pion yield. Solid curve: the final yield of free pions; dashed curve: the number 
of real quanta in the initial state. 

The effect of the enhancement on the pion yield is finally shown in the right panel 
of fig. P for all five events, using the same procedure as in the quench scenario. Based 
on these calculations one would expect a rather broad enhancement of the pion yield 
at low momenta momenta, p < 250 MeV/c, corresponding to kinetic energies below 
about 150 MeV. This enhancement is fairly modest, amounting to perhaps 50%, but 
might be visible as an excess relative to a thermal extrapolation based on the yield at 



23 



higher energies. In addition, there appears a fairly robust peak near the upper edge 
of the enhanced region. Its magnitude is probably difficult to predict reliably, since it 
is sensitive to the amplitude and persistence of the induced oscillations in the chiral 
order parameter. Furthermore, although the peak location can be pin-pointed rather 
well within the model calculation, it should be kept in mind that the determining 
quantity is the a mass which is not well established. In fact, the a field is introduced 
as merely a simplistic representation of the observed rather broad enhancement of 
the elementary isoscalar cross section. Consequently, it might be too naive to expect 
well localized peak to appear in the observed spectra. 



7 Concluding remarks 

In the present work, we have discussed the quantum field treatment of free particles 
exposed to an arbitrary time- dependent external field that may become supercritical 
so that spontaneous pair creation occurs. Though of general nature, this problem 
is particularly relevant for current efforts to explore the possibility for forming so- 
called disoriented chiral condensates in high-energy nuclear collisions, because the key 
effect is expected to be the amplification of soft pionic modes by either a transient 
supercritical field or the eventual oscillatory relaxation of the order parameter near 
its vacuum value. 

For the treatment of the time- dependent quantum-field problem, we adapt a 
method introduced by Combescure [0; this method was also employed in Ref. 



14]. It has some resemblence with the the treatment of the corresponding quantum- 



mechanical problem (where the goal is to obtain the wave function for a forced os- 



cillator) introduced by Husimi We have developed the general formalism and 
illustrated it both for instructive idealized cases and cases of direct DCC relevance. 
By divorcing the specific initial state from the behavior of the evolution operator, 
the method makes it easy to treat the problem for any given time dependence of 
the frequency, iJ^it). In particular, the expected final number of quanta in a mode, 
i^/, is simply given in terms of the initial occupancy z/q by means of an amplification 
coefficient that depends on iJ^if) but not on z/q. In addition to being convenient, 
this form is also instructive, as it brings out the fact that the vacuum fiuctuations 
and thermal fiuctuations contribute in a democratic fashion, both being subject to 
the same amplification by the factor X^. In particular, the result makes it easy to see 
how spontaneous particle production is possible even when starting in the vacuum 
state. 

As a supplement to the quantum-field developments, we also showed how the 
same key results can be obtained entirely within classical mechanics by performing a 
suitable analysis of how the results depend on the initial phase of the state. Although 
we suspect that this novel method is in fact mathematically equivalent to the quantal 
treatment, we have only demonstrated the correspondence by means of numerical 
calculations. The fact that processes that are inherently quantal in nature, such as 
spontaneous production in supercritical fields, can be treated accurately within the 
confines of classical physics may prove to be practically useful, since it is significantly 



24 



easier to treat the field dynamics classically. 

Moreover, the formal discussion also makes it possible to understand the quantita- 
tive limitations of a semi-classical treatment of the field equation. This is important, 
since by far most of the dynamical calculations in the context of disoriented chiral 
condensates have in fact been carried out classically. The developed classical method 
for obtaining the amplification coefficients may possibly be generalized and brought to 
bear on the evolution of the entire field in order to extract more reliable information 
from essentially classical calculations. 

Finally, employing time-dependent effective masses obtained from dynamical cal- 
culations within the linear a model, we have illustrated the degree of amplification 
that can occur in possible collision scenarios. In particular, a significant enhancement 
of certain modes is expected to occur as a consequence of the quasi-regular behavior of 
the chiral condensate, as it relaxes after its initial agitation far away from its vacuum 
value. In a rather extreme scenario, where the cooling rate emulates a scaling expan- 
sion in three dimensions and a transient quench is generated, the soft pion modes are 
enhanced by well over an order of magnitude. A more cautious scenario, where nei- 
ther the initial chiral restoration nor the subsequent expansion is as extreme, yields a 
modest but persistent enhancement of the soft modes as well as a larger preferential 
enhancement at a finite pion energy related to the mass of the effective a meson. This 
latter effect is a parametric resonance phenomenon whose appearance in the context 
of disoriented chiral condensates has already been recognized by several groups (see, 
for example, Refs. [|, [23[|). 

The presented applications serve primarily as illustrations only and we have lim- 
ited the discussion to the enhancement of the one-particle spectra. It has been sug- 
gested that the time dependence of the effective mass may produce a substantial 



signal in the two-particle correlations and recent studies have explored the pos- 



sibility of using this effect as a DCC signal [^. The present treatment may provide 
a helpful tool for further explorations of this interesting prospect. 

Once the time dependence of the effective mass has been calculated, in a suitable 
model, the present formalism makes it possible to calculate the resulting effect on the 
pion dynamics with due account taken of the complications arising from the inherent 
quantum nature of the system. In this regard, we expect that the present work will 
be of some practical utility. In particular, it might be possible to use it as a means for 
making a perturbative estimate of the quantum enhancement. Of course, it would be 
desirable to extend the treatment to the more general case where the effective mass is 
endowed with a spatial variation as well; this problem presents an interesting future 
challenge. 

Helpful and stimulating discussions with Andrei Krzywicki are gratefully acknowl- 
edged. This work was supported by the Director, Office of Energy Research, Office of 
High Energy and Nuclear Physics, Nuclear Physics Division of the U.S. Department 
of Energy under Contract No. DE-AC03-76SF00098. 
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Appendix: Coherent states 



We here recall some of the features of coherent states that may be especially useful 
in the context of the present study. 

One would ordinarily like to associate a free field with a specific coherent state, 

1;^) = eEjx^"^l-Xk^K.]\Q) = J^^eEnXnO^l\0) , (86) 

where = exp(— | J2k IXkH- Thus |x) is normalized to unity, 

ixlx) = KK' eS->^-^- . (87) 

Furthermore, a^\x) = Xk,\x) and = xIx'k.' ixlx') ■ For free motion the 

time evolution of |x) is governed by dtXn = —i^nXn- We also note that 

1 r2~ 

{xlQnlx) = ^={x\a^ + al\x) = W — Rex^ , (88) 



{x\Ph.\x) = y(^I"1*~"'cI^) = \/'2^Klrax^, (89) 



so that Xk = \J^K/'2qn + Vv2^Pk- Thus, if we write Xn = \^exp{-ir]K), we have 



/ LJ cos 1]^ , = -J2u^uj^smT]^ . (90) 



We also note that when uj"^ is constant in time, then the coherent state ( ^6]) is a 
solution to the equation of motion, idt\t) = H\t), where \t) = exp(X)|0) with 

X{t) = xa^ - X*a = i{pq - qp) ■ (91) 



We have here used the relations (l§-|y) to introduce the canonical operators q and p 
(which are independent of time). The new parameters (which are real) are related to 
the original (complex) parameter x{t) by 

X{t) = q{t) + ^ Pit) = {a) (92) 

and they represent the expectation values of the canonical operators, q{t) = {t\q\t) 
and pit) = {t\p\t), as the notation would suggest. We then note that the evolution of 
the state is given by \t) = expi—iHt)\to) when H is time- independent. It possible to 
show that the resulting state is still of coherent form and that the parameters q and 
p satisfy Hamilton's equations, 

q = p , p = —uj'^q . (93) 

Thus, in summary, when the Hamiltonian is of harmonic form, H = p'^/2 + uj'^q'^/2, 
then the coherent state \t) = exp(z(pg — qp)\0) is a solution to the equation of motion 
provided the parameters satisfy the corresponding Hamilton equations, q = p and 
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p = —uj'^q. This result does not depend on the sign of but rehes only on its 
constancy. 

The above coherent state ( ]86| ) is written in terms of the operators associated 
with the diagonal basis. The state can equally well be written in terms of the original 
operators Ok which annihilate quanta having a definite momentum k, 

\x) = e^'*t'^'*°'^~^^°'^l|0) = e^k[Xk4-Xk«kl|o) . (94) 

The parameters Xk are easily obtained. With «, > we find, 

Xk=/c = ' ^k=-K = :^[^'^ + ■ (95) 

A coherent state of effective (or 'squeezed') quanta (which is what would ordinarily 
be used as an initial state), is also a coherent state in terms of the original free 
representation, 

1^) ^ e*^'-**^|0) = e^'^'-^-'^lO) ^ \x) , (96) 
where Xk denote the coefficients pertaining to the free representation, 

, = 4±2S; = (^)*Rex + ,(^)*imi, (97) 

and we have used that — X* P = X^'^ ^ X^cx- Moreover, it has been assumed that 
the effective frequency is positive, > 0, as will be hte case if the quasiparticle modes 
are in thermal equilibrium. This relationship can be rewritten on a simple democratic 
form, 

X. = h^^ + ^p^ ^ x. = Vy^'^ + tis^^- (9^) 



where and are the expextation values of the canonical operators in the diagonal 
representation, as given by the (real) coefficients in the trigonometric expansions 

f|T2lO). 
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